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Abstract 

We introduce a flow-dependent version of the quadratic Steiner tree problem in the plane. 
An instance of the problem on a set of embedded sources and a sink asks for a directed tree 
T spanning these nodes and a bounded number of Steiner points, such that /(e) |e|^ is 

e6-B(T) 

a minimum, where /(e) is the flow on edge e. The edges are uncapacitated and the flows are 
determined additively, i.e., the flow on an edge leaving a node u will be the sum of the flows 
on all edges entering u. Our motivation for studying this problem is its utility as a model for 
relay augmentation of wireless sensor networks. In these scenarios one seeks to optimise power 
consumption - which is predominantly due to communication and, in free space, is proportional 
to the square of transmission distance - in the network by introducing additional relays. We 
prove several geometric and combinatorial results on the structure of optimal and locally op- 
timal solution-trees (under various strategies for bounding the number of Steiner points) and 
describe a geometric linear-time algorithm for constructing such trees with known topologies. 

Keywords: Power-p Steiner trees. Network flows. Mass-point geometry, Wireless sensor networks 



1 Introduction 

Given a set of points Z in a normed plane (M^, |1 • 1|) and a real number p > 0, the geometric power-p 
Steiner tree problem (or geometric p-STP) seeks a finite set of points 5' C (tlie Steiner points) 
and a tree T = {V{T), EiT)) = (Z U 5, E{T)) such tfiat ^ \\u ~ v\f is a minimum. For p > 1 

uv£E(T) 

tfie input to p-STP must include a strategy for bounding the number of Steiner points, without 
which a minimum solution may not exist. When p — 1 we obtain the classical Steiner tree problem 
[51 111), which has been extensively studied under the rectilinear and the Euclidean norms. Soukop 
|13) was the first to explore the notion of non-linear networks from a topological point of view; he 
realised its importance as a model for the design of transportation or communication systems. The 
operations research community has studied a very similar problem in the form of the non-linear 
multi-facility location problem; see for instance [5]. The p-STP, in the form given above with || • || 
the Euclidean or rectilinear norm, was introduced by Ganley and Salowe in [SJ (5 for its application 
to VSLI routing algorithms. 

The p = 2 case, which we refer to as the quadratic STP, is particularly important for transportation 
problems [17] and for some wireless network problems. In the latter case this is because most of the 
energy of the network is utilised during data transmission, and, furthermore, energy consumption is 

*This research was supported by an ARC Discovery Grant. 
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proportional to the transmission distance raised to an exponent a £ [2, 5]; see for instance [H]. Since 
a, the so-called path loss exponent, is greater than 1, adding a relay between any two communicating 
nodes will lead to a reduction in total energy consumed. This leads us to an effective method of 
reducing power consumption in wireless sensor networks through relay augmentation. In free-space 
a can be shown to be exactly 2 [8]; moreover, the a = 2 case is important in some real- world 
sensor network applications when constructive interference applies, such as in beamforming and 
communication through corridors [T^. We therefore solely address the quadratic STP in this paper, 
but with an additional flow component that makes for a more realistic model of relay augmentation 
in wireless networks. 

In order to make the above discussion more rigorous we note that the energy consumed by nodes 
in a wireless sensor network when a single packet of data is transmitted over a distance r in free 
space is /3r" -I- pi-oc = f3r'^ + Pica where /3 is a constant and p^cc is the energy required to receive a 
packet. Since Pioc is usually small relative to we simplify and normalise the energy consumption 
to E(itw) = jutip, where u is the transmitting node, v is the receiving node, | • | is Euclidean distance, 
and uv denotes the edge (u,v). In order to model power consumption we must include the rate 
of data flow between u and v, say f{uv), to get P{uv) — f{uv)\uv\'^. We assume that there is an 
additive flow function on the nodes, so that the difference between the flow rate entering a given 
node and leaving the node equates to the supply rate at the node. It is therefore a simple matter to 
calculate f{uv) for any m, v once we have the supply rates at all sensors (sources) and the topology 
of the network has been given. 

The central problem of this paper is referred to as the flow- dependent quadratic Steiner tree problem 
(FQSTP). An instance of the FQSTP has an input of n -I- 1 points in the plane - i.e., n sources and 
one sink - and a supply rate at each source. We are also given a strategy for bounding the number 
of Steiner points. The output is a set of Steiner points (satisfying the given bound) in the plane 
and a tree T interconnecting all nodes such that the sum of P(m, v) over all u, v is minimised. We 
can view the FQSTP as a weighted version of the quadratic Steiner tree problem, where the weights 
are the flows on the edges. A related flow-dependent problem is the Gilbert arborescence problem 
(GAP) [71 [Tl] , which also has an additive flow function at the nodes but where the cost of an edge 
uv is w{fuv)\uv\, where w is some increasing concave function. Similar to the FQSTP, the GAP has 
applications in telecommunications and transportation networks. The computational complexity of 
the FQSTP is unknown; in fact, this is true for the complexity of the geometric power-p Steiner tree 
problem for any p besides 1. Ganley [3] observes that the methods used for proving NP-hardness 
of the classical geometric Steiner tree problem fail for p > 1 because of the lack of the triangle 
inequality; however, when the number of Steiner points is bounded by some constant strictly less 
than |Z| — 2 then the quadratic Steiner tree problem can be shown to be NP-hard by reduction from 
the geometric dominating set problem. Berman and Zelikovsky [1 show that the graph version of the 
p-STP (where the Steiner points are restricted to being vertices of a given graph) is MaxSNP-hard. 

Our main contributions in this paper are an analysis of some of the important geometric properties 
of optimal (or locally optimal) solutions to the FQSTP, and a new linear time geometric algorithm 
for the construction of locally minimal trees. This algorithm can be used as a component in an 
exact solution to the FQSTP, and we also believe that it will lead to effective pruning techniques 
for such an exact solution (similar to those used in the program GEOSTEINER for solving the 
Euclidean and Rectilinear STP [H]). We provide a formal definition of the FQSTP in Section 
[2l In Section |3] we describe various strategies for bounding the number of Steiner points in order 
to ensure that a solution exists. Then, in Section |4l we state and prove a number of structural 
and geometric results on locally minimal solutions to the FQSTP for each of the afore-mentioned 
bounding methods. Finally, Section [5] presents our geometric linear-time algorithm for constructing 
locally minimal trees of a given topology. 
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2 Problem definition and notation 



Let Z U {2;bs} be a set of n sources and a sink (or 6ase station) zbs embedded in the plane; we 
assume that Z ^ $ and z-qs ^ Z. We assume that there is a supply w{zi) > associated with 
each member oi Zi & Z and a demand w{zi) associated with ^bs- These supplies and demand 

ztez 

determine the flow on the resulting tree. Let S be any finite set of additional points in the plane, 
referred to as the Steiner points. A directed tree T spanning Z U S U {zbs} is a flow-dependent 
quadratic Steiner tree (FQST) if and only if 

1. Every edge e of T is labelled by a positive real number /(e), its flow, 

2. T is directed towards Zbs such that every node z of T, except the sink, has exactly one outgoing 
edge e+ {out-edge) with flow /+ = fi/'t) '^^i*^! possibly some incoming edges {in-edges) with 
sum of flows /~. The sink has no out-edges, but has at least one in-edge, 

3. For the sink we have /^g = and /^g = w{zi), 

4. For each source z we have ft — fz = w{z), 

5. For each Steiner point s, ft = , 

6. The cost of T is L{T) := ^ f{ei)\ei\^, where E{T) is the edge-set of T. 

ei^E(T) 

The objective of the FQSTP is to minimise L{T) over all FQSTs. The decision variables for this 

problem arc the number and locations of the Steiner points and the topology of the tree intercon- 
necting all points. An FQST minimising L{T) will exist if and only if l^j is bounded, and below 
we discuss various strategies for doing so. An optimal solution will be referred to as a minimum 
flow-dependent quadratic Steiner tree (MFQST). Throughout this paper we assume that w{zi) = 1 
for all Zi G Z, but all our results can be generalised to any positive w(zi). 

It should be noted that points (3)-(5) in the definition of an FQST define an additive flow function 
on the nodes. This means that the flows from the sources are, in some sense, independent of each 
other; this fact leads to a sharing of some properties between MFQSTs and shortest path trees. The 
wireless network analog of this is that data aggregation (for instance compression) does not take 
place at the nodes. 

For any node z the neighbours of z incident to the in-edges of z will be called z's in-neighbours, 
and we have a similar definition for the out-neighbour of z, which we sometimes refer to as the local 
sink of z. The degree of a node is the number of edges incident to that node, and its in-degree is 
the number of in-edges incident to it. We denote an edge or a line segment connecting points u and 
V by uv, and its Euclidean length by \uv\. The familiar notation \\u\\ is used for the length of u 

considered as a vector, i.e., ||u|| = ^u^-\-u'^ where u = {ux,Uy). 

Any tree network T interconnecting some or all of the nodes of Z[JSU{zbs} induces a tree topology 
T, which is simply the labelled graph corresponding to the network T. If, in T, every z € ZU {zbs} 
is of degree one and every s G 5 is of degree larger than one then T is a full topology. Since Steiner 
points are never of degree one in a MFQST, the edge set of any T induced by an MFQST on 
Z U S U {^Bs} can be partitioned such that every member of the partition induces a full topology. 
A tree is called degenerate if and only if it has at least one edge of zero length. 
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3 Bounding the number of Steiner points 



As stated before, an MFQST on a set of nodes will exist if and only if there is a bound on \S\. To 
see this suppose first that there is no bound on IS*]. We can reduce the cost of any FQST on the 
given nodes by adding one or more degree-two Steiner point to any edge; see Fig. [1] Note that 

ItAT^I \ ^ I ZiZ^ I f(wU^ 

the total cost of such a (straight-line) path is (p + 1) ( — IjTy ) fi''^'^) — 1^ — where u,v are 

the end-points of the path, and the Si are equally spaced degree-two Steiner points (in Fig. [1] and 
throughout, the sources and sink are shown as filled circles and Steiner points are open circles). 
Since, if there is no bound on \S\, we can keep adding degree-two Steiner points, there is no optimal 
solution; in fact, the cost of the tree will tend to zero as the number of Steiner points increases. 
On the other hand, if \S\ is bounded then there are a finite number of possible tree topologies 
interconnecting the nodes, with each topology obtaining a unique minimum (as we show in the next 
section). Therefore an optimal solution must exist. 




Figure 1: Total cost of path is 



p+1 



A popular method of bounding the number of Steiner points in the power-p Steiner tree problem, 
p > 1 , is to restrict the degrees of Steiner points to be greater than two [H |3l HI [13] . It is clear 
that we then have an implicit upper bound ofA:=|S'|<n — 1. Alternatively we could introduce 
a cost on the Steiner points. This latter bound is the one that is most relevant to the application 
motivating this paper, namely wireless sensor network deployment [151 1181 [T^ . In this scenario the 
Steiner points correspond to transmitting relays, which have an (often significant) associated cost. 
In practice, there may also be a need to impose edge-length bounds on the network, but we do not 
study such bounds in this paper. 

Let T be an FQST. In summary, we consider the following three |S'|-bounding strategies. 

1. The degree bound stipulates a fixed value (j) > 3 such that (p < degs for any Steiner point s in 
T. 

2. The explicit bound stipulates a fixed upper-bound of k for \S\. It should be clear then that 
jS*! = /c since adding a Steiner point always leads to an improvement in cost. 

3. The node-weighted version of the problem assigns a weight c > to every Steiner point. The 
objective of the node-weighted FQSTP is then to minimise Lc{T) :— /(e) |ep -f els']. 

eeE{T) 

4 Properties of locally minimal FQSTs 

Here we consider combinatorial and geometric structural properties of MFQSTs and locally minimal 
(with respect to a given topology) FQSTs. In the first subsection we describe general properties 
that hold true under any |S'|-bounding strategy. We then look at properties relating specifically to 
degree-bounded FQSTs, and finally we look at properties for FQSTs that are not degree-bounded. 
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4.1 General properties 



Suppose that we are given a set of points Z U {zbs} embedded in R^, a set S — {si, s^} of (free) 
Steiner points and a topology T interconnecting Z U S* U {^bs}- For any u — {ui, ...,Ufc) G M^'"' 
(where each Ui G R^) let Tu be the tree with topology T that is obtained by embedding Steiner 
point Si at location Ui for every i. We wish to minimise L{T^) = /(e)|ep over all u € R^*^. 

Such a minimum is unique, when Tu is non-degenerate, by the strict convexity of L(Tu) (noting that 
it is a sum of strictly convex functions). A locally minimal tree with respect to 7" is a tree Tu of 
topology T minimising L{T-^). Clearly any MFQST is locally minimal with respect to its topology. 

Now let Z = {zi, Zn) be any set of points in the Euclidean plane, and suppose that we associate 
a mass fi with Zi for each i G / := {1, . . . , n}. We use the familiar mass-point geometry notation 
fiZi to refer to this weighted point, and we denote the centre of mass of the system of points 
M = {fiZi i G /} by C{M) or (by a slight misuse of notation) C{fiZi, fnZn)- 

Ml 



^ O fn+l Z„+i 



fn 

Figure 2: s = C(/iZi, /„+iz„+i), where /„+i = ^fi 



fiZi 




Proposition 1 Suppose that a Steiner point s in a locally minimal FQST T has in-neighbours 
zi,...,z„ providing respective flows and out-neighbour Zn+i such that fn+i ■— /(ej") = 

fi (see Fig. l^. Then s ^ C{fiZi, fn+iZn+i)- 

i<n 

Proof. Let / = {!,..., n}. Suppose that we perturb the Steiner point t units away from s in the 
direction of the unit vector u and let the resultant tree be Tq. Then 

ijjuit) := L{To) = ^fi\\s + tu~ Zi\f + /„+i||.s + tu- Zn+i\f 
attains a minimum at t = 0. But 

■ip'uit) = 2^ S + tu- Zi) + 2/„+i.(u, S + tu- Zn+l) 

and therefore 

V'L(O) = 2^/,.(u,s-z,) + 2/„+i.(u,.s-z„+i) =0. 
iei 

Therefore 



, X! - ) = fn+l{Zn+l - s)) . 



U 

iei 
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Since this holds for any unit vector u we get 



Therefore 



^ft{s- Zi) ^ fn+l{Zn+l - S). 



f i + /n+l 

and the resuh foUows since this is an expression for the weighted mean. ■ 

According to the method of mass-point geometry [2j, any point C(/izi, fnZn) can be constructed 
geometrically by recursively merging masses and subdividing line segments into appropriate ratios. 
Consider for instance the construction of C{fiZi, f2Z2, (/i + 12)23) where /i,/2 > and the Zi are 
any three points in the plane, as illustrated in Fig. [31 We first merge fiZi and f2Z2 into the point 

\z\Z\ 2 1 f\ 

(/i + f2)zi2 where zi 2 is a point on the line segment Z1Z2 such that — - = - — . Merging 

\Z2Z1^2\ h 

(/i + f2)zi.2 and (/i + f2)z3 yields the point s — — since zi^2 and 23 have equal masses. 

We extend this merging method in Section [5] in order to construct locally minimal FQSTs for more 
general topologies. 



(/l+/2)2l,2o^"" 



/l 

Figure 3: Merging masses to construct a Stciner point s. Here t = — —\ziZ2\ 

J I + h 



The reasoning behind the following corollary should now be obvious. 

Corollary 2 Any Steiner point in a locally minimal FQST lies at the mid-point of its out-neighbour 
and the centre of mass of its in-neighbour s, where masses are assigned to the neighbours of the Steiner 
point as in Proposition^^ 

Fig. m shows an example of a locally minimal tree on four sources. Note the length ratios that 
result from the previous corollary; in particular, t — C{lz, 3s2), and si lies at the midpoint of t and 
z'. 

Next we show that the converse of Proposition [T] is also true. The proof is similar to Ganley's proof 
for quadratic Steiner trees [3]. 

Theorem 3 An FQST is locally minimal if and only if every Steiner point lies at the centre of mass 
of its neighbours. 
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Figure 4: A locally minimal tree on four sources and a sink 



Proof. We only need to prove sufficiency. Let T be an FQST with Steiner points si, Sp, p > 0. 
Let the in- neighbours of s.; be x\,..., x'l\ providing flows //,■••, ff\ and let the out-neighbour of Si 
be x'j^ which receives a flow of /{efj = //• Furthermore, suppose that every Steiner point of T 

3<di 

is at the centre of mass of its neighbours, so that 

Si = ; — 

where all three sums are for 1 < j < d^. Therefore 



Let A be the p x p square matrix containing the coefficients of the Si in the previous equation; 
note that some (or all) of the xj or x'^ may be Steiner points. This creates a system of linear 
equations As = b, where s = (si, Sp), and b is a constant derived from the fixed locations of the 
sources and the sink. By Proposition [T] one of the solutions to this system must produce a locally 
minimal tree. However, by Equation ([1]) the diagonal entry in the ith row of A is 2^^//, which 
has magnitude strictly larger than any other entry of row i. Hence A is diagonally dominant and 
therefore non-singular by the Levy-Desplanques theorem; see for example [lOj. ■ 

The next lemma is an interesting result on adjacent Steiner points in locally minimal FQSTs. The 
lemma is also useful for proofs in Sections 14.21 and [5] Let sq, si be any two adjacent Steiner points 
of a locally minimal FQST T, where si,zi, are the neighbours of so; so,Zj+i, ...,Zp+i are the 
neighbours of si; and Zp+i is the local sink. Once again let fi be the flow associated with the edge 
incident to Zi. Let / — {1, ...,p+ 1}, J = {1, for j < p, and J — I — J. Let Cj — C{{fiZi}i^j) 
and Fj = /j , and similarly for C-j and F-j. 
ieJ 
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Zi 



Figure 5: The points Cj, so, si, C-j are coUinear 



Lemma 4 With the above notation and definitions, the points Cj , Sq, Si,C-j are collinear and sub- 
divide the line segment CjC-j into a ratio Fj : Fj : Fj. 

Proof. The lemma is illustrated in Fig. [S] Note that si = C {FjSq, F-jCj) and sq = C{FjCj, F-jSi), 
from which the result follows. ■ 

4.2 Properties of degree-bounded MFQSTs 

We next examine some properties of degree bounded MFQSTs, where the given lower bound on 
degree is > 3. Let x d Z U S he a node in a locally minimal FQST T, with in- neighbours zi, ...^ Zp 
and local sink Zp+i. Let / = {1, ...,p + 1} and let J C I\{p + 1}. A J-split (or just split if the 
context is clear) of x introduces a Steiner point s' such that the neighbours of s' are x and every 
Zi, i G J, and the neighbours of x are s' and Zi for every i £ /\J. We assume that, after splitting, 
s' and all other Steiner points are relocated to their optimal positions relative to the new topology. 
Let the resultant tree be denoted by Tq. A J-split is beneficial if L{To) < L{T). Most splits are 
beneficial, but not all; see Fig. [6l 

In this section we use the existence of beneficial J-splits to show that imposing a lower bound </> > 3 
on the degree of Steiner points implies there is also an upper bound of 2^ — 3 on the degree. We 
need the following definition and lemma before we can show that beneficial splits of any size can 
always be found. 

We define two edges in a network to be overlapping if they are both incident with a common node 
and every point of one edge lies in the other. 

Lemma 5 Let T be a non- degenerate FQST with a pair of overlapping edges, where either T is not 
degree-bounded, or T is degree bounded and the common node of the two overlapping edges is not a 
Steiner point of degree (j). Then T is not an MFQST. 
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Figure 6: A {1, 2}-split of s is not beneficial since C(lzi, IZ2) — C(lz3, 824) — s 



Proof. Suppose that T is as specified but is also an MFQST. Suppose that the two overlapping 
edges are ei — uiv and 62 = U2V with \uiv\ > \u2v\. Then, regardless of the directions of flow, we 
can replace edge ei by U1U2 and thereby reduce the total cost of T. This is a contradiction. Note 
that the edge replacement is allowed in the degree-bounded case because the only node-degree that 
decreases is that of v, which, by assumption, is not a Steiner point of degree 0. ■ 

In the next proposition we let ds denote the in-degree of s. 

Proposition 6 Let T be a non-degenerate locally minimal FQST containing a Steiner point s such 
that ds > <j) — 1 if T is degree-hounded (with degree hound (p), and ds > 1 otherwise. Furthermore, 
assume that no pair of overlapping edges in T have s as their common node. Then there exists a 
heneficial J-split of s for any 1 < | J| < ds. 

Proof. Suppose that the in- neighbours of s are zi,...,Zd^ and let J C {!,..., dg}. We argue by 
induction on the cardinality of J. Clearly the result holds for | J| = 1 since C{fiZi) = Zi ^ s for 
any z < ds by non-degeneracy; also, the result holds for |J| — ds since C{fiZi, fd^ZdJ / s by 
Corollary [21 Next suppose that \J\ — d for some 1 < d < ds — 2 such that the J-split of s is 
beneficial, and let a, a' G {1, ds}\J be distinct. Let Jq — J U {a} and let Ji = J U {a'}. Note 
that if Cj„ = Cjj — s then the three points Za, Za' and s are coUinear, with Za or Za' lying between 
the other two (again by Corollary [2|) . This contradicts the assumption that s has no incident pair 
of overlapping edges. Therefore, either the Jo-split or Ji-split of s is beneficial and is of cardinality 
d+1. m 

We now have the following four results for the degree bounded problem. The two corollaries also 
hold for Euclidean quadratic Steiner trees; see [13]. 

Proposition 7 An MFQST T with degree hound (p has <f) < deg s < 2(j) — 3 for every Steiner point 
s. 



Proof. Suppose, contrary to the proposition, that T contains a Steiner point s with deg s > 20 — 3. 
Note that s has no incident pair of overlapping edges, by Lemma[5] Any J-split of s with | J| =0—1 
will produce two Steiner points each of degree at least <j). By Proposition [5] at least one choice of 
such J must be beneficial, which contradicts the minimality of T. ■ 



Corollary 8 When = 3 every Steiner point will be of degree exactly 3. 
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Proposition 9 Every source z in an MFQST T with degree bound 4> is of degree at most 0—1. 
Moreover, if z has degree equal to (j)—! and if C denotes the centre of mass of z and its in-neighbour s, 
then z lies at the midpoint of C and the out-neighbour of z. 

Proof. The reasoning is similar to the previous proposition, except that when using J-splits there 
is no lower bound on the in-degree of z. ■ 

Corollary 10 Every source in an MFQST with (f) = Z is of degree at most two. Moreover, if a 
source has degree equal to 2 then it is collinear with its neighbours (see for instance Fig. 



4.3 Properties of MFQSTs that are not degree-bounded 

In this subsection we consider some properties of the node-weighted and the explicitly bounded 
versions of the FQSTP. 



Proposition 11 The angle between any in-edge and the out-edge of a given node in a node-weighted 
or explicitly bounded MFQST is at least 90° (see Fig. 

Proof. If we assume that 9 < 90° as in Fig. [5] then path z^, Zj is shorter (quadratically) than path 
Zi, X, Zj (by Pythagoras' theorem). Therefore, since there are no restrictions on the degree of nodes, 
we can replace edge z^x with edge ZiZj, resulting in a tree with lower cost. ■ 

In fact, by repeatedly swapping edges we can ensure that every in-edge-out-edge angle is strictly 
greater than 90°, although this may cause the sink to obtain a large degree 




Figure 7: An MFQST with (f> — 3 with a source of degree 2. 




Figure 8: The angle between an in-edge and an out-edge must be obtuse 
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Figure 9: Nodes in an MFQST can acquire large degree 



Proposition 12 Steiner points of a node-weighted or explicitly bounded MFQST can achieve a 
degree of \Z\ + 1. 

Proof. For the explicitly bounded case, suppose that sources are located on a circle as in Fig. [9l 
with a single Steiner point. For this FQST to be locally minimal, the Steiner point must be located 
near the centre of the circle, but slightly towards the sink. Clearly, by Proposition 111! this tree is 
also globally minimal - in particular, if any of the sources has degree greater than 1 in the MFQST 
we get a contradiction to the angle condition. Therefore the Steiner point for the MFQST with 
explicit bound fc = 1 in this instance has degree \Z\ + 1. 

For the node-weighted case, let \Z\ = n, and suppose the sources and sink are located on a unit 
diameter circle such that all sources lie within an arbitrarily small neighbourhood of the antipodal 
point to the sink. Then for any number of Steiner points the minimum cost of the network will be 
close to that for a network with a single source of weight n. Let the cost of each Steiner point be 
n/2 — e for a small e > 0. Then it is clear that the MFQST has exactly one Steiner point, located 
close to the centre of the unit diameter circle. By the same argument as in the previous case, this 
Steiner point has degree n + 1. ■ 

Next we present a few properties relating to the node-weighted version of the FQSTP. Recall that 
this version utilises a modified cost function, namely Lc{T) — /(e) |ep + c\S\ where c > is 

eeE{T) 

the cost of a Steiner point. The first two results are necessary conditions for an FQST to be an 
MFQST, and are based on the number of degree two Steiner points that lie on any given straight 
line path of an optimal tree. 

Proposition 13 Letp be the number of degree-two Steiner points located on a path P with endpoints 

f luv^ I uv I 

u,v of a node-weighted MQFST. Then p{p + 1) < < (p + 2)(j' + 1), where c is the cost 

c 

of a Steiner point. 



Proof. Note that the nodes on P are collinear and equally spaced along the segment uv. The 

-, „. f(uv)\uvP , . , f(uv)\uv\^ f(uv)\uv\'^ 

cost of path P is — — - — - + cp. We therefore need — — - — - + cp < — — - — - + c{p + 1) and 

p+l p+1 p + 2 

\- cp < h c{p — 1), from which the result follows from simple algebra. ■ 

p+l p 



Corollary 14 In a node-weighted MFQST every edge e satisfies |e| < 




where c is the cost of 
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a Steiner point. 



Since we are not directly bounding k in the node-weighted version, it would be helpful to deter- 
mine an upper bound B for k in terms of c. This would immediately lead to an exact algorithm 
for calculating node-weighted MFQSTs: the first step would be to iterate through all topologies 
interconnecting the given sources and at most B Steiner points. A locally optimal solution for every 
topology could then be calculated using either the algebraic or geometric methods described the 
next section, and the cheapest tree selected. The complexity of this method, however, would be 
prohibitive for large problems (i.e., large n or small c) and effective methods of pruning the number 
of viable topologies would be needed in these cases. We first prove the following result in order to 
bound the "length component" of the cost of a node-weighted MFQST. 

Lemma 15 Let T he an explicitly hounded MFQST on n sources, with bound k. Then L(T) > 
(n + fc + l)-i^|z,ZBs|'. 

i<~n 

Proof. Consider a single unit of flow from source Zi. A cheapest path from Zi to zbs must have 
a cost of at least that of the path P, where P contains all Steiner points and remaining sources 
arranged collinearly and equally spaced between and zbs- In this case the cost of P would be 

— ■ — . ■ 

n + k + 1 

Let T be a minimum spanning tree on ZU{zbs}, and suppose that we add the minimum number of 
degree-two Steiner points to T such that no edge is longer than ^J^, where / is the flow on the edge. 
Note that this can be done greedily, and therefore the construction is of polynomial complexity. Let 
the resultant tree be denoted by BST{Z), i.e., the beaded spanning tree on Z. 

Lemma 16 Let ZU{zbs} be any set of sources and a sink, and let Topt be a node-weighted MFQST 
on these nodes. Then Lc{Topt) < Lc{BST{Z)) . 

In the previous lemma it may be possible to construct a spanning tree (in polynomial time) that 
provides a tighter bound than a minimum spanning tree. An improvement of this bound, or the 
bound in Lemma 1151 would most likely lead to a better value of B in the next proposition. 

Proposition 17 Suppose that a node-weighted MFQST T contains k Steiner points. Then k < 
— ^Lc{BST{Z)) — (n + fc + 1)^^ l-^i^BSpl? which can be rewritten (by solving a quadratic equa- 
tion in k, in order to make k the subject) in the form k < B where B does not contain k. 

Proof. The result follows directly from the previous two lemmas after noting that k ~ — {Lc{T) — L{T)). 
_ c 



5 A Geometric Linear-time algorithm for fixed topologies 

In this section we describe a geometric linear time algorithm for constructing an MFQST with a 
given full topology, but where all Steiner points are assumed to have degree 3. This is a key step 
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in any general algorithm for constructing MFQSTs over all possible topologies. It can potentially 
be combined with an exhaustive search, along with appropriate pruning methods, to build an exact 
method for finding MFQSTs (along the lines of the method underlying GEOSTEINER [16]), or can 
be combined with appropriate heuristic search techniques as part of an approximation algorithm. 
As mentioned before, for many variants of the problem the restriction of the degree of Steiner points 
(especially to degree 3) is a very natural one. An algebraic linear time algorithm does exist in the 
form of a solution to the system of diagonally dominant linear equations discussed in the proof of 
Theorem [3] (see also [3] for the related algorithm without flow), but this algorithm reveals very little 
directly about the structure of locally optimal FQSTs. 

The general strategy of the algorithm is similar to Melzak's algorithm for constructing fixed topol- 
ogy Euclidean Steiner trees; see for instance [TT]. We begin by finding two sources that are adjacent 
to a single Steiner node in the given topology, and replace the pair of sources by a single quasi- 
source whose location and mass can be explicitly computed. This procedure is repeated recursively: 
at each stage there exists a Steiner point adjacent to two nodes, each of which is either a source 
or quasi-source; these nodes along with the Steiner point are replaced by a new quasi-source. The 
procedure continues until there are no sources left, and only a single quasi-source in the tree. The 
position and mass of this quasi-source allows one to construct the Steiner point adjacent to the sink, 
and then a backtracking procedure allows one to construct each of the remaining Steiner points in 
turn. 

Notation: As before, let Z ~ {zi, . . . , z„} be the set of sources for a full MFQST T, and let zbs be 
the sink. We think of each source and the sink as being a vector in representing the position of 
the node in Cartesian coordinates. Associated with each source Zi is a mass u'(zi), representing the 
amount of flow from that source. As in the earlier sections, we assume that w{zi) = 1 Vi G {1, . . . , n} 
under a suitable choice of units, however the algorithm can easily be adapted to situations where 
different sources have different flows. Let S = {si, . . . , s„_i} be the set of Steiner points of T. We 
associate a mass with each Steiner point Si and the sink additively; for example, the mass of each 
Si is the sum of the masses of the two nodes whose out-edges are the in-edges of Si. Nodes that are 
not Steiner points (including the sink) will be referred to as terminals. 

We now define the concept of a quasi-source. Given an MFQST containing a Steiner point adjacent 
to two terminals (neither of which is a sink) , a quasi-source is a new terminal that replaces these 
three nodes, such that the remaining Steiner points of the resulting MFQST on this reduced set of 
terminals are in the same locations as in the original tree. More formally, let T be a full MFQST 
with terminals Z U {^^bs} (where zbs is the sink), Steiner points S = {si, . . . ,s„^i} each with a 
given mass, and full topology T. Suppose we are given zi,Z2 G Z and si £ S such that zi and 
Z2 are each adjacent to si, and let v be the third node (other than zi and Z2) adjacent to si. 
This is illustrated in Fig. [TUl Let Ti be the subtree of T in which zi, Z2 and their two incident 

T T, T' 

Figure 10: The transition from T to T' (via Ti). 

edges have been removed, and let 7i be the topology of Ti. Note that if we treat si as a terminal, 
then 7i is a full topology. Given a point qi with mass w{qi) let T' be the MFQST with terminals 




13 



{qi, Z3, . . . , Zn} U zbs, topology 7i (where node Si has been relabelled qi), and where the weights 
of the Steiner points in T' are the same as their weights in T. Then qi is said to be a quasi-source 
replacing 2:1,22 and si, if the Steiner points of T' are in the same locations as the corresponding 
Steiner points of T and si — C(w{qi)qi,w{v)v). The transition from T to T' (via Ti) is illustrated 
in Fig. [TOl 

It is important to note that the above definition involves a slight abuse of notation - strictly 
speaking, T' is not an MFQST since the mass assigned to a quasi-source no longer has an obvious 
interpretation in terms of fiow. It is convenient, however, to treat T' as though it is an MFQST. 
Because the weights of Steiner points do not change when we introduce a quasi-source, it follows 
that such 'MFQSTs' that include quasi-sources as terminals do not necessarily have the additive 
property for all Steiner points (but, by definition, still satisfy the centre of mass properties). 

In the lemmas that follow we show that for a given T we can always find a quasi-source q, such 
that the position and mass of q depend only on the known masses of nodes in the tree and the 
locations of the two terminals it is replacing. We distinguish three methods of constructing such a 
quasi-source, depending on the nature of the two terminals it replaces. 



Lemma 18 Let zi and Z2 be two sources ofT, both adjacent to a single Steiner point s. Let q be 
the midpoint of zi and 22 (ie, q — {zi + Z2)/'2) and let w(q) = 2. Then q is a quasi-source replacing 
z\ and Z2. 



Proof. Note that q is the centre of mass of zi and Z2, and w{q) = w{zi) ^ w{z2)- Let v be the third 
neighbour of s, other than zi and 22- It follows that the centre of mass of zi, 22 and v is the same 
as the centre of mass of q and v. Hence the result follows. ■ 

Lemma 19 Let q and z be two terminals of T , a quasi-source and source respectively, both adjacent 
to a single Steiner point 82- Let si be the Steiner point (not in T) adjacent to the two terminals 
replaced by q in a previous minimum quadratic flow-dependent Steiner tree Ti. Let wa = w{q) and 
wi = w{si) (in Ti). Define the point q\ as follows: 

qi:=q + ■ ■ (2 - q) 

Wo + Wi + WqWi 

with mass 

Wo-\-Wi-\- WqWi 

w(qi) := ■ . 

Wq -f Wi 

Then qi is a quasi-source replacing q, 2 and 82- 



Proof. The aim is to show that the point qi as defined in the statement of the lemma is a suitable 
quasi-source for replacing q, z and 82- Consider the construction illustrated in Fig. 1111 Let v be the 
third node adjacent to 82 other than si and 2. Note that v is either a Steiner point or the sink, and 

w{v) — Wi + 1. 

In Ti, si is the centre of mass of q and 82- Hence q, Si and 82 are collinear. Let x be the point 
where the line through q, si and sq intersects the line segment vz. By mass-point geometry, 

M = ^—. (2) 

\xz\ wi + 1 

Furthermore, we can consider x to have an associated mass of w{x) = wi + 2. It follows that 
-r^^il = l£i. and li-SsI = which, in terms of ratios, gives: 

\qsi\ : IS1S2I : |s22:| = wi{wi -\- 2) : wa{wi -f 2) : wiWq. (3) 
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Figure 11: Construction for the proof of Lemma [T9l The weights of some points are shown in 
parentheses. The single and double arrows are used to indicate parallel lines. 



Now, let L be the line through q parallel to vz, and let t be the point where the line through si and 
V intersects L. Let a be the point where sit intersects qz. Since Atqsi ^ Avxsi it follows from ^ 
that 1^ = "'!("' 1+^) and hence, by (12]), = 75 — /"^ , Since, Atqa ^ Avza, we now obtain: 

\xv\ 2ii;o(^f 1+1) ' ' ^ ' — \zv\ 2wq{w-i + 1) ' ^ ' 

M = . (4) 

\qz\ wi + 2'Wo{wi + 1) 

Let b be the point on qz such that S2b || sia. Since Aqasi ^ Ag6s2, we obtain from ^ and (|3]): 

\ab\ Wo 



\qz\ wi + 2wo(wi + 1) 
Similarly, let c be the point on qz such that S2C || xz. Since Acs2q Azxq, we deduce that: 

\cz\ WiWq 

\qz\ {wi+wo){wi+2)+wiWo' 



(5) 



(6) 



Equations (H)), ^ and ©, give us the locations of points a, b and c respectively. Now let qi be 
the intersection of the line through v and S2 with qz. Since Abcs2 ^ Aazv, we obtain: 

_ \ ab\ + |6gi 



|6c| 

Equation ([T]) allows us to compute |6(7i|, resulting in the location of qi as given in the statement of 

\q1S2 _ \bqi 
\qiv\ ~ \a.qi\ 



the lemma. Finally, we note that AbqiS2 ^ Aaqiv, hence = |^2li from which we obtain the 



mass of qi as stated in the lemma. 



Lemma 20 Let qi and 52 be two quasi-sources of T , both adjacent to a single Steiner point S3 . 
For i — 1,2, let Si be the Steiner point adjacent to the two terminals replaced by qi in a previous 
minimum quadratic flow- dependent Steiner tree Ti, and let woi ~ w{qi) and Wi — w{si) (in Ti). 
Define the point 53 as follows: 

W02W2(W1+W2) , . 

13 qi H 7 : 7- 7 ; r 92 - qi) 
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with mass 

WiW2{woi + W02) + W0lW02{wi + W2) 

wiqa) ■■= 7 — - — w — X • 

(Wi +Woi)[W2 +W02) 

Then (73 is a quasi-source for T, replacing qi, (72 cind S3. 

Proof. The proof is similar to the proof of Lemma [T^ Consider the construction given in Fig. [T^ 
where again v is the third node adjacent to S3 other than si and S2, and w(v) = wi + W2- let Li be 




Figure 12: Construction for the proof of Lemma BUI 

the hne through qi parallel to VS2^ and let L2 be the line through q2 parallel to wsi; for i = 1, 2, let 
ti be the point where the line through Si and v intersects Li . As in the proof of Lemma 1191 can 
compute the points a and 6, where qiq2 intersects siii and 52^2 respectively, and the points c and d 
on qiq2 such that S3C and s^d are parallel to va and vb respectively. This again allows us to locate 
the point q^ at the intersection of the line through v and S3 with (7192, with location and mass as 
given in the statement of the lemma. ■ 

We now describe the algorithm for locating the Steiner points of T by successively replacing pairs 
of terminals by quasi-sources using the above three lemmas. 



Algorithm MFQST 

Input: A set Z oi n sources, a sink zbs, and a full topology T for Z U {zbs}- 

Output: A minimum flow-dependent quadratic Steiner tree T for ZU{zbs} with topol- 
ogy T, along with its cost L{T). 

1. For each Zi & Z set w{zi) = 1; for each Steiner point Sj of 7" compute w{sj), via 
additivity; set u'(zbs) = n. 

2. Find a Steiner point adjacent to two terminals (each of which is a source or quasi- 
source) and replace the three nodes by a new quasi-source, using Lemmas [THl [TH] 
and [201 Update the Steiner tree. Repeat until all Steiner points have been replaced. 

3. The Steiner tree will now contain only two terminals, a quasi-source and the sink. 
Use recursive back-tracking to determine the positions of the Steiner points (in the 
reverse order to the order of replacement in the previous stage) where each Steiner 
point is at the centre of mass of its two neighbouring nodes. 
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4. Once all quasi-source replacements have been undone, the tree T with all of its 
Steiner points will have been correctly constructed. The cost of T is computed 
using node weights and edge lengths. 

Before proving the correctness of Algorithm MFQST, it is helpful to illustrate the running of the 
algorithm with a simple example, shown in Fig. 1131 where all locations are given in Cartesian 
coordinates. The example contains four terminals: three sources zi — (0, 0), Z2 = (2, 4), Z3 ~ (11, 5), 




Figure 13: Example of the construction of an MFQST using the algorithm. 

and a sink zbs = (H, !)• The full topology of the tree is shown in unbroken lines. In Step 1, the 
weights of the two Steiner points are computed by additivity: u;(si) = 2 and 'w{s2) = 3. In Step 2 
we replace each Steiner point and two adjacent terminals by a quasi-source. The first quasi-source 
qi replaces si and the two sources zi and Z2. Hence, by Lemma |18[ qi = (zi -I- Z2)/2 = (1,2) and 
w{qi) = 2. In the resulting tree, the remaining Steiner point S2 is adjacent to qi and Z3. We replace 
these three nodes by a new quasi-source q2 where, by Lemma [T9l q2 — qi + (zs — 9i)/2 — (6,3.5) 
and w{q2) = 2. This concludes Step 2. For Step 3 we determine the positions of the Steiner points 
in reverse order to Step 2. The Steiner point S2 lies at the centre of mass of 92 and zbs, hence 
S2 = 1^2 + f ^BS = (9, 2). Similarly, si lies at the centre of mass of qi and S2, where, for this branch 
the relevant mass of S2 corresponds to the flow on the edge S1S2 in the final tree (ie, 2). Hence 
si = i(7i + = (5, 2). This completes Step 3, giving the MFQST T with L(T) = 102. 

Theorem 21 Given a set of terminals Z and a corresponding full topology T , Algorithm MFQST 
correctly computes a MFQST T for Z. Furthermore, the algorithm runs in time 0{n). 

Proof. We claim that, in Step 2 of the algorithm, if the current tree contains Steiner points then 
there is at least one Steiner point adjacent to two terminals, neither of which is the sink. If the tree 
contains only one Steiner point then the statement is trivial, while if the tree contains more than 
one Steiner point, then this follows from the observation that the subtree induced by the Steiner 
points contains at least two vertices of degree 1. Hence we know that such a Steiner point exists at 
each stage, and the correctness of the algorithm easily follows from Lemmas [TSl flQl and [20l 

For the running time, note that the order of replacing Steiner points by quasi-sources can be pre- 
determined by running a depth-first search on the topology, which can be done in linear time. The 
construction of each quasi-source can be done in constant time, using the formulas in Lemmas 1181 
1191 and 1201 hence Step 2 of the algorithm runs in linear time. Similarly, Step 3, determining the 
position of each Steiner point, requires only linear time. ■ 

Finally we note that the methods in this section can easily be extended to allow Steiner points of 
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both degree 2 and 3. The inchision of higher degree Stcincr points, however, appears to significantly 
increase the complexity of the problem, and may require a different approach. 



6 Conclusion 

In this paper we introduced a flow-dependent version of the quadratic Steiner tree problem in 
order to model optimal relay deployment in wireless networks, specifically networks that transmit 
in free space or in situations where constructive interference applies. We described some structural 
geometric properties of locally minimal solutions to the problem, including properties relating to 
the degrees and locations of Steiner points. We did this under various strategies for bounding the 
number of Steiner points. Finally, we described a new geometric algorithm for constructing locally 
minimal solutions. The algorithm is based on the mass-merging method of mass-point geometry, 
and runs in linear time, matching the fastest known algebraic algorithm for the problem. 
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